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ABSTRACT 

In the ultraviolet (UV), Type la supernovae (SNe la) show a much larger diver- 
sity in their properties than in the optical. Using a stationary Monte-Carlo radiative 
transfer code, a grid of spectra at maximum light was created varying bolometric lu- 
minosity and the amount of metals in the outer layers of the SN ejecta. This model 
grid is then compared to a sample of high-rcdshift SNe la in order to test whether 
the observed diversities can be explained by luminosity and metallicity changes alone. 
The dispersion in broadband UV flux and colours at approximately constant optical 
spectrum can be readily matched by the model grid. In particular, the UVl-b colour 
is found to be a good tracer of metal content of the outer ejecta, which may in turn 
reflect on the metallicity of the SN progenitor. The models are less successful in re- 
producing other observed trends, such as the wavelengths of key UV features, which 
are dominated by reverse fluorescence photons from the optical, or intermediate band 
photometric indices. This can be explained in terms of the greater sensitivity of these 
detailed observables to modest changes in the relative abundances. Specifically, no sin- 
gle element is responsible for the observed trends. Due to their complex origin, these 
trends do not appear to be good indicators of either luminosity or metallicity. 
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1 INTRODUCTION 

Type la supernovae (SNe la) are some of the most impor- 
tant tools for current cosmological studies. Following the 
disco very that the ir peak magnitudes could be standard- 
ised (|PhiIlipsil 1993ft . thei r use enabled the discovery that the 
unive rse is accelerating |Riess et al.lll998l ; IPerlmutter et all 
Il999l h More recent supernova observations, combined with 
other constraints from the cosmic microwave background 
and baryon acoustic oscillations have established that we 
live in a flat unive r se with a matter content of Q,m ~ 0-27 
(jAstier et alJ 120061 : IWood-Vasev et alJ 120071 : iKessler et all 
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|2009| ; ISullivan et alJliuTll ) and the remaining 73% made up 
of dark energy, the nature of which is currently unknown. 

SNe la studies have measured that w, the dark energy 
equation of state parameter is consistent with w = — 1 to 
6.5% with some studies at high-z even beginning to place 
constraints on whether w = w (z) (|Riess et alj |2007h. In the 
future, we will have to observe at higher redshifts in or- 
der to find supernovae from younger times in the universe 
to improve on these dynamical measurements. As such our 
surveys will either have to switch to the IR or probe the 
rest-frame UV, a region of the spectrum that has been less 
extensively explored in the local population. 

The UV spectra of SNe la have long been thought to 
prob e the region where metallicity effects would be impor- 
tant (|Hofhch et all 1 19981 ; iLentz et al1l200ut ) due to the vast 
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number of metal line transitions in this region. Many of the 
photons in this region are absorbed, mostly by iron group 
elements, in what is referred to as metal line-blanketing; 
however, the effect that metallicity has on the level of th e 
continuum flux in this region is debated (|Sauer et al.ll2008h . 

One example of the effect that progenitor metallicity 
may have is that in a higher metallicity progenitor, the pro- 
duction of neutron-rich iso topes such as 54 Fe a nd 58 Ni is 
favoured compared to 56 Ni (|lwamoto et al.|[l999h . This will 
be reflected not only in the spectra, but also to some degree 
in the broad-band light curves. A SN la in the local universe, 
where metallicity is high, will on average have a lower lumi- 
nosity than a high-redshift SN la due to the different 56 Ni 
content and hence a d ifferent light curve stretch , as shown in 
iTimmes etal] |2003l 'l; lHowell et all (|2007l . l2009l ). The abun- 
dance ratio of stable iron-group elements to radioactive 56 Ni 
has been proposed as an addition al parameter for the stan- 
dardi sation of SN la light curves (|Mazzali fc Pod siadlowski 
l2006h . 

Until recently, our understanding of the role of progen- 
itor metallicity in SNela data was limit ed by the paucity 
of observed UV spectra. Recognising this, lEllis et al.l l|200gl . 
hereafter, E08) secured high quality Keck spectra for 36 in- 
termediate redshift (z ~0.5) SNela at maximum light drawn 
from the Supernova Legacy Survey (SNLS); th ese o ptical 
spectra appropriately probe the rest-frame UV. E08 noted 
a significant diversity in their UV spectra which could not be 
attributed to dust. Importantly, they found the variations in 
their UV data, as characterised by colours derived directly 
from their rest-frame spectra, did not correlate with the 
light-curve stretch. They also showed that the wavelengths 
of specific UV features showed phase-dependent shifts. The 
dispersion in their UV colours was claimed to be larger than 
could be accounted for meta l-dependent models available at 
the time (|Lentz et al.ll200(ih thus opening the possibility of 
an additional explanation for the diversity. 

Recently, Hubble Space Telescope (HST) and Swift ob- 
servatio ns have begun to e xplore the rest-frame UV of local 
SNe la. iFolev et al] (|2008l l used archival HST and Interna- 
tional Ultraviolet Explorer (IUE) data to show that a par- 
ticular ratio of UV flux correlates strongly with absolute V- 
band magnitude for 6 objects with spectra near maximum 
light: brighter supernovae have lower valu es of the ratio. This 
is a different result from that claimed in lEOSl which saw no 
correlation with supernova brightness. 

In a more rece nt study us i ng 21 intermediate redshift 
SNe la ~0.25), IFolev et al] (|2012T ) repeat the UV flux 
ratio analysis and find a different relation between absolute 
V-band magnitude and the ratio value. In this case, brighter 
supernovae still show lower values of the ratio value, but the 
slope of the relation is very different. IFolev et al] (|2008l ) ex- 
plored the use of the UV ratio as a luminosity indicator for 
light curve standardisation with some degree of success. The 
potential link between UV properties and intrinsic luminos- 
ity would have important implications for future cosmolog- 
ical studies. 

Swift data ha ve been used by iBrown et all (|2010h and 
iMilne et all (|2010h to obtain an overview of the spectral be- 
haviour in the UV. In the near-UV filters (2600-3300A and 
3000-4000A), the normal sub-class of SNe la shows a high 
degree of homogeneity, while the subluminous and the pecu- 
liar SN 2002cx-like groups show large differences. Absolute 



magnitudes at maximum brightness are correlated with the 
optical decay rate and show a scatter similar in size to that 
obtained with optical data. However, in the mid-UV (2000- 
2400A) the scatter is much larger (~ 1 mag) , indicating pos- 
sible metallicity-driven effects in this part of the spectrum. 

Recent HST o bservations of 12 Hu bble-flow SNe la at 
maximum light by ICooke et al] (|201ll) show that the dis- 
persion from a mean spectrum increases as wavelength de- 
creases, and is largest in the UV region of the spectrum. 
They attribute this to the larger number of metal absorp- 
tion lines in the UV compared to the optical. The same 
increas e in d ispersion is seen at higher redshift to the same 
degree (|E08h so they conclude that this must be an intrinsic 
feature of the supernova and not due to evolutionary effects. 
A larger study of UV spectra at maximum from HST is un- 
derway (Maguire et al. 2012, in press). A large degree of 
diversity is also seen in the UV photome t ry an d spectra of 
four supernovae discussed in I Wang et al] (|2012h . The paper 
concludes that more detailed modelling of supernovae in the 
UV is required. 

Optical studies have shown that supern ova properties 
depen d on the properties of the host galaxy. lHamuv et al] 
|2000t ) first showed that the higher mass galaxies preferen- 
tially host dimmer supernovae compared to brighter super- 
novae which were associated with younger stellar popula- 
tion s in late-type g a laxies . 

ISullivan et aTJ (|2010l ) examined SNe la subdividing 
their sample by host properties. They found that in more 
massive galaxies, or in those with a lower specific star- 
formation rate, the SNe la were on average ~ 0.08 mag 
brighter than that of SNe la in other gala xies after cor- 
rectio n for light-curve stretch and colour. ISullivan et al] 
l|2010l ) suggested that the difference they observe may be 
due to the metallicity of the host galaxy as more mas- 
sive galaxies tend to be more m etal-rich; however thi s ap- 
pears at odds with the results o f ITimmes et al] (|2003T l and 
iMazzali fc Podsiadlowski (|2006h . 

This study is thus motivated by the need to reconcile 
the conflicting deductions regarding the observed diversity 
in the UV spectra derived from earlier work We exploit a 
wide range of models parameterised by both bolometric lu- 
minosity and metallicity to see if we can explain the obser- 
vations with these two variables alone. The model dataset is 
presented in Section [2. II and the optical data sample which 
we use for comparison is described in Section ^. 21 In Section 
[3]we compare measurements of various UV properties of the 
model and data samples. Our results are then discussed in 
Section [3] and a summary of our conclusions presented in 
Section [5] 

2 DATASETS 

2.1 Model Dataset 

We calculate model spectra in order to study the influence 
of metallicity and luminosity on the UV spectra of "normal" 
SNe la around _B-band maximum. The respective r a diativ e 
transfer models are based on those in ISauer et al. | (|200Sf ). 
A two-parameter grid of models is set up by changing the 
bolometric luminosity of the models and the metallicity in 
the outer ejecta to study the effect of these parameters on 
the UV. 
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2.1.1 Code and Model Input 

To calculate the synthetic spectra, we ran a stationary 
Monte-Carlo radiative transfer code, which has successfully 
been used to model phot ospheric spectra of numerous SNe 
la in spherica l symmetry dMazzali fc Lucvl l 19931 ; lLucvl[i"99i 



iMazza li 2000: Stchl e et a.1.1 120051 ) . including rare UV obser- 



vations (|Sauer et alj 20081 ) 



The version of the code used here simulates a SN atmo- 
sphere above a lower boundary ("photosphere"), and takes 
as input data the location of the photosphere, r p h or v p h', the 
time passed from explosion onset, t; an abundance stratifica- 
tion; and a density profile. Here, we use the d ensity profile of 
the S N la delayed detonation model WDD3 ()lwamoto et al.l 
1999). The profile is scaled by the code, assuming homolo- 
gous expansion to time t, i. e. for each mass element r = vt 
(where r is the distance from the centre of mass and v is 
the velocity imparted at explosion). Thus either r or v can 
be used as spatial coordinates. Apart from the parameters 
mentioned, the code allows the user to set a bolometric lu- 
minosity, Lt i, for the final output spectrum. 

From the photosphere, which is located at an adjustable 
w p h, energy packets of continuous black body radiation 
(Jj" = B v T p h) are assumed to be emitted into the atmo- 
sphere. The simulated photon packets undergo Thomson 
scattering and interactions with atomic absorption/emission 
lines, which are treated in the Sobolev approximation. A 
downward branching scheme ensures a good approximation 
to the bound-bound emissivity. Radiative equilibrium is en- 
forced by t he construction of the Monte-Carlo simulation 
(|Lucvll999l ). Bound-free processes are not simulated, as lines 
togethe r with Thomson scattering dominate the opacity in 
SNe la ijPauldrach et al.lll99rj ; fSauer et al.ll2006h . 

The excitation and ionisation state of the matter are 
calculated from the radiatio n field statistics using a modi- 
fied n ebular approximation ()Mazzali fc Lucvlll993l : iMazzalil 
l2000h . In this approximation, the gas state in each radial 
grid cell is mostly determined from a radiation tempera- 
ture Tr and a dilution factor W. Tr corresponds to the 
mean frequency of the radiation field and W parametrises 
its energy density (for given Tr). We iterate the state of 
the gas and the radiation field in turn until the Tr values 
within the atmosphere are converged to the per cent level. 
Within these iterations, T p h is automatically adjusted so as 
to match the given luminosity Lboi- This adjustment com- 
pensates for the reabsorption of radiation which occurs when 
packets re-enter the photosphere via back-scattering. After 
the iterations have converged, the final spectrum is calcu- 
lated by formal integration of the transfer equation (|Lucvl 
1 19991 ). using a source function derived from the Monte Carlo 
statistics. 



2.1.2 Model Grid 

All synthetic spectra in this paper are derived from a model 
developed for SN 2005c f, for which UV/optical data near 
maximum is available ()Bufano et al.l 120091 : iGaravini et al.l 
120071 : 1 Wang et all 120091). SN2005cf wa s a bright supernova 
(M B = -19.39 JPastorello et alj|2007t)) with a light-curve 
width Amis = 1.10 (jPastorello et all 120071 : IWang et all 
120091 ) which correspo nds to a light- c urve stretch of s = 0.93 
using the relation in IConlev et al.l l|200Sf ). Since SN2005cf 
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Figure 1. Synthetic spectrum from our ejecta model (red line) 
compared to an observed UV- optical spectrum (black line) taken 
0.9 da ys before B maximum teufano et alj|2009l ; Ig aravini et al.l 
l2007f ). 



was a bright SN and produced a large amount of 56 Ni 
(« 0.7JV/©) we chose to use as the density structure of model 
WDD3 from llwamoto et all |l999t) . The bolometric luminos- 
ity of this model is \og(Li, i / Lq) — 9.6. 

We have used an abundance tomography approach to 
model the UVOIR spectrum of SN2005cf 0.9 days before 
maximum. The method of successively constraining the 
abundances in deeper and deeper layers usin g a temporal 
sequen ce of spectra was first introduced by IStehle et al.l 
(|2005l ). In order to constrain properly the highest veloc- 
ity material in the outermost ejecta before trying to model 
a spectrum at maximum, we started with a spectrum ob- 
tained —7.8 days before B-band maximum. The maximum 
brightness spectrum has a phase of —0.9 days relative to B- 
band maximum. We assume a B-band rise-time of 19.5 days 
JConlev et alj|2006h and use t = 18.6 days. In order to op- 
timise the fits to the data, two additional zones were intro- 
duced on top of those at maxim um and pre-maximum ve- 
locities, as in lSauer et al.l (|2008l ). This is not unexpected as 
strong high-velocity f eatures have been noted in SN 2005cf 
(|Garavini et al.l [20071 ). We optimised the match to the UV 
spectrum (A < 4000 A) while maintaining the best possible 
fit in the optical. 

This model is reproduced in Figure [T] A summary of 
the chemical composition in each zone is given in Table [T] 
where IME stands for intermediate mass elements (those 
with atomic number 9 - 20) and IGE for iron-group elements 
(atomic number greater than 20). 

From Table [1] we see that the outer two layers are dom- 
inated by unburnt carbon and oxygen with very little IME 
and almost no IGE. In the pre-maximum layer, there is still 
some C/O material remaining and slightly more IGE mate- 
rial, but the shell is dominated by material that has been 
burnt to IME. The shell at the photospheric velocity at max- 
imum is still dominated by IME, but now a significant frac- 
tion of the shell is made up of IGE. 

We used this technique to create models for supernovae 
with different bolometric luminosities. The range of lumi- 
nosities was chosen to reflect the spread of observed SNIa 
properties. We also used different underlying density profiles 
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Table 1. A summary of the major element groups in each shell 
of the model for SN 2005cf. X is the mass fraction of C/O, inter- 
mediate mass elements (IME) and iron group elements (IGE). 



Zone 


"ph 


"pre- max 


%ut,l 


^out^ 


Velocity 


10750 


13100 


16000 


19500 


X(C/0) 


0.3 


4.0 


70 


92 


X(IME) 


63 


92 


18 


8.2 


X(IGE) 


37 


3.9 


3.2 


0.2 



Table 2. An overview of the model parameters for the models 
log( L ho ;/L w ) = 9.2 9.4 9.6, 9.75 with r? = 1. The W7 model is 
from lNomoto et al.l lll984t) and the WDD1 and WDD3 models arc 
from jlwamoto et al.lll999h . The given velocities are in km s . 



9.2 9.4 9.6 9.75 



Density Structure W7 WDD1 WDD3 WDD3 

v ph 6783 8539 10750 11250 

upre-max 8366 10406 13100 13710 
«out,l 16000 
«out,2 19500 



to reflect the different energies and Ni production at the 
bolometric luminosities \og(Lboi/Lo) = 9.2,9.4,9.6,9.75. 
The velocities of the and tipre-max shells move depend- 
ing on luminosity (Table 2), resulting in the masses of indi- 
vidual elements being scaled for the whole model. However, 
within each shell the relative abundances of the SN 2005cf 
model are preserved as described. In order to ensure that 
we produce realistic models, we compare the model output 
spectrum to observed SNe la to ensure we are producing 
spectra which match the continuum levels in the UV while 
maintaining normal optical spectra i. e. are not members of 
the over-luminous of under-luminous subclasses. This pro- 
cess is summarised in Tableland the spectra are displayed 
in Figure 

We then used these 4 luminosity bins to create a set 
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Figure 2. The r\ = 1 spectra for \og{L bo i/ Lq) = 
9.2, 9.4, 9.6 , 9.75 . Also marked are the three top-hat filters in- 
troduced in E08, UV1, UV2 and b (dashed lines) and the Bessel 
U filter (solid black line). 
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Figure 3. The r) = 5,2,1,0.5,0.2,0.1,0.05 spectra for 
log(Lj, ;/ £q) = 9.6. The features Ai and A2 are also marked. 

of models with varying metal contents. We do this by scal- 
ing the metallicity in the pre-maximum and two outer shells 
with respect to the best-fitting model: the photospheric shell 
remains unchanged. To generate the sequences, we multi- 
plied the abundances of all the elements with atomic number 
Z > 20, i. e. heavier than calcium, by a factor 77, which is 
allowed to take the values 77 = 0.05, 0.1, 0.2, 0.5, 1, 2, 5. The 
mass fraction X of element E thus becomes X(E) = riX(E)o 
where X(E)o is the mass fraction of the element in the 
SN 2005cf best-fit model at the expense of unburnt C/O. 
This provides us with a grid of 4 x 7 models for our analysis, 
but we have excluded the model where log(i(, ; / Lq) — 9.75 
and 77 = 0.05 as the optical spectrum did not look nor- 
mal, leaving 27 models for analysis. A metallicity sequence 
is shown in Figure [3] for \og(L bo i/ L i3 ) = 9.6. For all values 
of 77, the optical spectra appear normal with relatively little 
dispersion. However, in the UV the dispersion between the 
models increases dramatically and diversity is also seen in 
the shapes and positions of features. This reflects the fact 
that metal line-blanketing effect is stronger in the UV. 

One caveat with our models is that in the red and in- 
frared the fits to data are less good. This is because of the 
crude assumption of a blackbody at the photosphere. The 
flux inside an SN la is non-thermal even in the inner layers 
fsee ISauer et al.ll2006l ). Flux redistribution within the in- 
ner parts of the simulated atmosphere leads to a sufficiently 
accurate radiation field in the atmosphere in the ultravio- 
let and blue regions of the spectrum, but in the red and 
infrared some flux excess usually remains in the synthetic 
spectra with respect to observations. This can be seen in 
Figure [1] where the model flux is higher than the observed 
flux from ~ 6300 A onwards. This means that our estimates 
of LboI may be somewhat larger than the real value when 
the red and IR are overestimated. Therefore, in order to 
compare models and data we extract Lb from both. 

2.2 Observational Data 

The observa tiona l data for this study are taken from those 
presented in |E08| . These SNe were discovered as part of the 
Supernova Legacy Survey (SNLS, ISullivan et al.l l201ll ). a 
real-time Type la supernova search based at the Canada- 



France-Hawaii Telescope (CFHT) which used g' r' i' z' -band 
observations to identify high-z SNe la. For more de- 
tails on the r eal-tim e SNLS target selection pipe-line see 
iPerrett et all |2010r ). 

[EOS! observed a selection of SNLS supern ovae on the 
Keck I Telescope, using LRIS ()Qke et al.lll995h to obtain a 
high signal-to-noise ratio in the rest-frame UV. Using host 
photometry obtained from the deep stack CFHT images 
in the u*</r'i'z'-bands, a best-fitting template galaxy spec- 
trum was used to remove contamination from host g alaxy 
lig ht in the SN spectru m. For more details on this see IE08i . 
or IWalker et al.l (|201ll ) which gives a detailed explanation 
of the application of this method to SNLS data obtained at 
the Gemini and VLT telescopes. 

It is imp ortan t to note that while the sample of super- 
novae used in |E08| as a whole was representative of the su- 
pernova population, the sub-sample of these spectra we use 
here are not because we apply a cut for rest-frame phase. In 
this study, in order to make a realistic comparison of mod- 
els to data, we included only supernovae with a rest-frame 
phase of ±3 days from B-band maximum. Additionally, we 
only considered "normal" SNe la with good lightcurve cov- 
erage so the B-band maximum magnitude could be calcu- 
lated, and with a spectrum reaching a minimum rest-frame 
wavelength of ^ 2700A. These cuts leave 9 objects. 

Within our sub-sample the mean light-curve stretch is 
s = 1.08, where 1.0 represents the fiducial "normal" SN la 
lightcurve. In fact, within our sub-sample, all but one of the 
objects have a stretch value > 1. The mean colour of our 
subsample is c = —0.027. T he light-curve fittin g was carried 
out using the SiFTO fitter (|Conlev et alj|2008f ). SiFTO also 
fits a B-band magnitude at maximum. To obtain Lb we 
convert magnitude to flux and the n to luminosity assu ming 
a flat cosmology with Qm = 0.269 (|Sullivan et al.|[201lT ) and 
Ho = 70kms- 1 Mpc~ 1 . 

In Figure [4] we plot the various models at 
\og(Lboi/L®) = 9.6, 9.75 and compare them to our 
observational data assuming all supernovae are placed at 
10 pc. The variations in metallicity appear to match the 
variations in the spectra of the IE08I sample. The spectrum 
lying below the models and the other observed SNe la with 
log(Lboi/ Lq) = 9.6 is 04Dlsk. This appears to have a low 
UV flux compared to the optical. As such, we expect this 
to be an outlier in some of our UV colour analysis (Section 
13.11) . None of the observed supernovae have luminosities 
which would correspond to the log(L(, ( / L@) = 9.2,9.4 
models; however SNe la are observed to have luminosities 
in this range (see Figure [TT] for example). 



3 ANALYSIS 

We now have a set of models and observed spectra which 
can be analysed in identical ways. As the B-band magni- 
tudes of the observed supernovae are well-known from their 
use in cosmology, we revert to using the B-band luminosity 
Lb instead of bolometric luminosity as this quantity is mea- 
sureable for both the models and observed supernovae. We 
can therefore make direct comparisons between the observed 
and model data for the UV diagnostics described in this sec- 
tion. It is possible to see from later plots that at constant 
Lboi, Lb does not vary linearly with 77. 
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Figure 4. For all values of r], the models with \og(Lb i/ Lq) = 
9.6 (top) and log(Lj io ;/L ( 7) ) = 9.75 (bottom) are compared to 
observed spectra from lEOSI . 

3.1 UV Colours 

We first examine the bro adban d UV colours at maximum 
light as first examined in |E08| . Their study used Bessel U 
filter as well as three top-hat filters they defined in the UV 
and optical (UV1, UV2 and a normalisation filter b; see 
Figure [2}. We do not perform any correction for supernova 
colour and remeasure the U, UV1, UV2, b magnitudes for 
the models and the observed spectra. The results are plotted 
in Figure [5] 

This shows that the broad trends observed within the 
lEOSl data at maximum are replicated in the models. The fil- 
ters U and UV2 do not show any strong trend with colour 
and B-band luminosity. The dispersion between models of 
different metallicity is small, and the trend between metal- 
licity and colour is not always linear. The UV1 filter shows 
that for the higher metallicity models the UV colour can be 
large for supernovae with lower luminosities. In general, the 
dispersion in the UV1 filter, which is the bluest of the three, 
is the largest. Given the smaller dispersion in the models and 
the large errors in the data, it is not possible to attribute 
the change in colour to metallicity, although this may be- 
come possible with more, better data. In the UV1 filter, the 
evidence that our less luminous SNe come from regions of 
higher metal content is stronger. 

As predicted above, 04Dlsk is the reddest object for 
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Figure 5. UV colours for the models and the data. The models are shown as the coloured squares and the data in the black filled circles. 
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Figure 6. The UV emission features Ai (left) and A2 (right) from lE08l plotted against B-band luminosity for both observed data and 
values measured from the models. The different colours denote the different metallicities in the model spectra. The observed data is 
divided into sub-samples based on the rest-frame phase of the spectra 



all three colours. It appears particularly red in the UVl-b 
bands implying that its ejecta are particularly metal-rich. 



3.2 UV Emission Feature Diagnostics 

Another property examined by |E08| is the position of two 
UV features described as Ai and A2 positioned at ~ 2900A 
and « 3100A respectively. The variation of these peaks with 
metallicity can clearly be seen in our models in Figure[3] The 
reason for this shift can simply be found in the fact that if 
the metal content is higher lines are effective at higher veloc- 
ities, and hence absorb at bluer wavelengths, progressively 
shifting the reemission peaks towards the blue. 

We use the same Gaussian-fitting method as described 
in lEOSl to measure the positions of Ai and A2 in the model 
data as well for the observed spectra, although we first ap- 
ply some smoothing to the observed data. As we only have 
model spectra at one phase (maximum) we choose to fur- 
ther sub-divide this sample due to the steep o bserv ed time- 
dependence of the position of these features jEOct ). These 
comparisons are shown in Figure [6] where instead of plotting 



Ai or A2 against phase, we plot them against Lb- The ob- 
served data are divided as —3 phase < 1 days (pre-max; 
open diamond); — 1 ^ phase < 1 days (max; filled circle); 
1 ^5 phase ^ 3 days (post-max; open circle). 

From the left plot in Figure[6]we see that at higher lumi- 
nosities there is a lower dispersion in the values of Ai in the 
models. The models show that Ai is roughly constant with 
luminosity in the lower metallicity models. We also see an 
approximately constant value of Ai within the observed data 
for all but 2 of the spectra with no strong phase-dependence. 

Figure [6] (right plot) shows that the position of the 
synthetic A2 peak emission varies strongly with both lumi- 
nosity and metallicity (see above). The observed data at 
±1 day agrees with the measurement in the models and the 
strong phase-dependence of the position of this feature is 
observed. At the highest luminosity, the position of the fea- 
ture is st rongl y blueshifted out of the range of measurements 
found bv lEOct This is because these features are strongly af- 
fected by the supernova velocities (see below) and at the 
highest luminosities the widths of Ai and A2 are very broad. 
This causes the two peaks to merge making the identifica- 
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4 DISCUSSION 



4.1 UV Colours 



Figure 7. The measured Ruv f° r both models and observed 
data. The different colours represent the different values of 7] for 
a given luminosity. 



tion of A2 difficult as the features are not Gaussian in nature 
(see Figures [2] and [3]) . 

3.3 UV Flux Ratio 

We can us e our data to measure the UV flux ratio of 
I Foley et al.l (|2008l ). Ruv, which is defined as 



Ri 



/a(2770A) 
/a(2900A) 



(1) 



where / A (2770A) and / A (2900A) are the median fluxes in 
bands of size ±20 A centred at these wavelengths. The re- 
sults are plotted in Figure [7] 

Figure[7]shows that for a given metallicity (77), the value 
of Ruv is approximately constant. The model trends agree 
broadly with the values of Ruv measured from the observa- 
tional data, but in this region of the spectrum the observed 
data is strongly affected by noise. We could only reproduce 
the correlation betw een Lb and Ruv a s would be expected 
from the results in iFolev et al.1 (|2008h if the brighter ob- 
served SNe la are in more metal-rich hosts. Instead the mod- 
els show that the diversity is driven almost entirely by the 
difference in metallicity. Some trend in seen in the data and 
the strength of this correlation is 2Ao using a Spearman's 
Rank correlation statistic assuming all points are equally 
weighted. We do not conclude anything from this because of 
the small sample size. 

We have also performed an analysis of Ruv after per- 
forming a colour- correction to th e supernovae using the 
SALT2 colour law (jGuv et al.ll2007f ). The change in the mea- 
sured value of the ratio with this correction in the order 
of 0.005 and so much less than the la error on the Ruv 
measurement based on the noise of the data. We have also 
examined what would happen shou ld the host subtraction 
of the galaxy as described in lEOSl be incorrect. We found 
that addition/subtraction of any galaxy template used in 
lEOSl with a u*-band flux of up to 20% that of the supernova 
would not change the Ruv measurement by more than the 
la error. As such we do not believe our Ruv measurements 
would be strongly affected by any inaccuracies in the host 
subtraction. 



We broadly replicate the wide range of colours observed in 
the UV with our models. We also show that dispersion in- 
creases towards the blue. This is due to the increasing extent 
of the metal line-blanketing which affects the whole of the 
UV region reflected in the opacity of the SN ejec t a increasing 
by 3 orders of magnitude from 4000A to 2000A l)Sauer et al.1 
|2006| ) . This increase in dispersion i s also consistent with the 
SWIFT data results iBrown et al.1 i|2010l ') 

The UV colours are the best diagnostic we have fea- 
tured for replicating properties seen by the whole observed 
data sample. This is because they measure over wide wave- 
length ranges and so are less affected by variations caused 
by the velocities in the models not matching the high-z data 
precisely. 



4.2 UV Peaks 

The two features designated as Ai and A2 in |E08| are not 
emission features in the traditional sense. The sheer number 
of species absorbing in these regions makes this extremely 
unlikely. The reason for the peaks in flux at ~ 2900 A and 
~ 3100 A could be due to two things, or most likely a combi- 
nation of both: a "window" in opacity at these wavelengths 
which allows more of the photons emitted from the photo- 
sphere to escape, or reverse fluorescence from species in the 
outer layers of the ejecta are generating these blue photons 
<|Mazzalill2O"O0l ). We can use the output of the model code to 
investigate what happens to photon packets emitted from 
the photosphere to try and differentiate between the two. 

Firstly, we look at what happens to packets emitted 
from the photosphere in a region of ±40 A around the wave- 
lengths of the measured values of Ai and A2. We find that 
at these wavelengths, for both features and at all luminosi- 
ties and metallicities, between 50% - 80% of the packets are 
reabsorbed back into the photosphere through backscatter- 
ing. Of the small number of packets left, they are mostly 
re-emitted at redder wavelengths. This means that an opac- 
ity window at low velocities can be excluded as the dominant 
effect. 

We can also look at photon packets which were origi- 
nally emitted in other regions of the spectrum, absorbed in 
some line transition and finally re-emitted within the Ai and 
A2 regions. This is illustrated in Figure [S] We see here that 
for Ai by far the most packets are gained from the red in 
reverse fluorescence processes, although this percentage de- 
creases as luminosity increases. The situation is similar for 
the A2 feature, although by the time we reach our highest 
luminosity model, the percentage is roughly 50% from the 
blue and 50% from the red. Reverse fluorescence was already 
mentioned as being prima rily responsible for the emerging 
UV flux by[M 

azza H (|2000l ): photons are fluoresced back to 
the UV at high velocities, and thus they can escape in low- 
oacity environments. 

If we also look at this measurement in terms of the num- 
ber of packets gained rather than the percentage of packets, 
we actually see that for both features the number gained 
from the blue increases with luminosity, which is to be ex- 
pected because of the increase in temperature. The number 
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Figure 8. The percentage of packets gained in the Ai (top) and A2 (bottom) regions which were originally emitted at redder (left) or 
bluer (right) wavelengths. 



of packets gained via reverse fluorescence, meanwhile, does 
not increase in the same way. For A2 the number increases 
in the lower two models and after that remains constant 
meaning that the lines in the optical have become saturated 
and this is limiting the number of photons re-emitted in this 
region. The number of packets for Ai does not saturate and 
so this is still affected by the abundances of the IME and 
IGE which absorb photons in the optical. 

We can also examine which ions cause the processes that 
shift flux into the Ai and A2 regions. We find that for both 
features, the predominant species which reverse fluoresce are 
IME, particularly Mgn, Sill and S 11, with some contribu- 
tions from the Fe-group elements, notably Cm. As these 
features are dependent on abundances of elements they will 
be affected by any differences between the velocity struc- 
ture and the abundance stratification between the models 
and the data. This can account for some of the discrepancy 
we see between the positions of the features. 



4.3 UV Ratio 

It has been established that a source of UV flux is a process 
of reverse fluorescence where red photons are absor bed and 
re-emitted at shorter wavelengths at high velocties (|Mazzalil 



l2000f ). This means that there have to be some metals in the 
material above the photosphere or there would be no UV 
flux at all. As the metal content is increased, the UV flux in- 
creases due to an increasing amount of reverse fluorescence; 
however this is in competition with an increased absorption 
of the UV photons by the metals themselves. The way this 
occurs in the regions of /(2770) and /(2900) determines how 
Ruv changes with metallicity at constant luminosity. 

We have plotted the values /(2770) and /(2900) in our 
models in Figure[9] In all cases, we see that the higher rj mod- 
els have less flux as there are more metals absorbing in these 
bands. In both bands we see a dispersion with metal content, 
which is larger for /(2900). This shows that the variation in 
Ruv is driven more by the variation in the /(2900) band. 

We used the MC code to identify the dominant ab- 
sorption lines in the regions probed by Ruv looking at ions 
with a large absolute change in absorption strength between 
rj = 0.05 and rj = 5 (Table [3]) or ions which have large mean 
absorption strengths averaged over the metallicity sequence 
(Table 13}. We can do this analysis for Ruv because of the 
small wavelength ranges for /(2770) and /(2900). No single 
element appears to be responsible for the change, although 
the /(2900) region appears to be dominated by variations in 
the Cm lines. We see that in the higher luminosity models, 
the dominant species change from singly-ionised to doubly- 
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Figure 9. Median flux values in the two regions probed by Ruv f° r our model (L,tj) sequences. 



Table 3. The species which show the largest absolute change in 
line strengths with metallicity in the two regions probed by Rjjv ■ 
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Table 4. The species which show the largest mean line strengths 
across all metallicitics in the two regions probed by Ruv ■ "~ " 
marks mean that the region contains a mix of singly- and doubly- 
ionised species of a variety of elements. 
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ionised reflecting the increase in temperature in the atmo- 
sphere. Figure [9] explicitly shows that effect of varying rj has 
on the two spectral bands. 

4.4 Focus on Individual Elements 

As well as creating sequences of models where the multiplier 
r\ affects all elements with Z > 20, we can use r\ to vary 
one element while keeping all others fixed. This allows us to 
examine whether this one element is the cause of any of the 
observed effects we see. Should any of this be caused by one 
element alone this would give us an important diagnostic to 
examine the physics of the supernova explosion itself. The 
species for which this analysis was carried out are chromium, 
stable iron, manganese, 56 Ni and titanium which are the 
elements with species with the strongest features in the near- 
UV. We summarise the effect of individual elements below: 

• UV Colours: Ti and Mn have no effect on the UV 
colours. Stable iron has little effect, except in the UVl-b 



colour where it shows strong non-linear trends with iron 
content at a given luminosity. The chromium and 56 Ni 
model sequences both show a large dispersion in the UVl-b 
colour as seen for all metals in Figure [5] above. 

• UV Features: The position of Ai is strongly depen- 
dent on the increased amount of Cr in the spectrum. This 
is expected from looking at which elements have caused the 
reverse fluorescence in this feature. The position of A2 is 
affected by the IGE produced in nucleosynthesis i. e. iron, 
nickel and chromium, but not manganese. We see that 
these fGE are small contributors to the reverse fluorescence 
flux which is dominated by IME, so Ai and A2 will be very 
sensitive to small changes in the relative abundances of all 
these elements. 

• UV Ratio: The sequences with varying stable iron and 
chromium content from nucleosynthesis all show large vari- 
ations at constant luminosity. This is expected as their ions 
were identified in Tables[3]and|3] There is diversity shown for 
all of the element sequences although to a lesser degree. This 
implies that this feature could vary significantly between su- 
pernovae and be dependent on velocity (which absorptions 
are shifted into the bands); temperature or luminosity af- 
fecting ionisation balanaces; and the abundances of the in- 
dividual elements relative to each other which may not be 
the same as for SN 2005cf. 

The plots detailing these results can be found in the Ap- 
pendix. 

4.5 Links to Other Observables 

4-5.1 Host Metallicity 

It is not possible to use our models to explain directly the 
observed trends with host properties. This is because our 
method of generating spectra does not separate out the effect 
of an increased metal content in the outer layers of the SN 
ejecta due to progenitor or environmental metallicity from 
any potential upmixing of elements which are burnt during 
the SN explosion. 

Hints on whether the environmental metallicity plays a 
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Figure 10. UV colours in the models and the observed spectra with the observed objects from lEoil divided into two sub-samples based 
on t heir host mass (top) and the host specific star formation rate (bottom). The values of the specific star formation rate (sSFR) are 
from lSullivan et al J 12013 ). The model sequences for different values of the metal content in the SN outer layers (rj) are again plotted in 
different colours. 



major role for the UV colours may directly be obtained from 
the observations by subdividing the observed sample accord- 
ing to metallicity indicators. In Figure fTOl we have done this, 
with the host galaxy mass and its specif i c star formation 
rate, sSFR, as measured in lSullivan et al.l ijil ) 1 n : as indica- 
tors. More massive galaxies with lower star-formation rates 
should generally be old, red ellipticals which have a high 
metallicity. In the upper panel of Figure 1101 there are two 
SNe in massive (high-metallicity) galaxies which are found 
to have law-Ruv / high-luminosity; however, there is also 
one SN in a high-mass host which shows an inverse tendency 
so we do not see any trend with host mass. Using the sSFR 
as a criterion (Figure [TOl lower panel), we see a similar re- 
sult. Fr 



4-5.2 Supernova Light-curve Standardisation 

iFolev et al.l (|2008l ) proposed that an observed correlation be- 
tween Ruv and the V'-band magnitude could be used for 
standardisation of supernova magnitudes. In Figure [TT1 we 
extend this into the B-bandand compare the supernovae in 
that sample to those from lE08l and our models, where the 
values of Mb from the models are generated via synthetic 



photometry. For the IFolev et alJ (120081) data, t he B-band 
magnitudes are taken from Altavilla et al.l |2004l ). The low- 
z data includes the measurement of Ruv nearest to max- 
imum for SN1980N, SN 1981B, SN1990N, SN 1991T and 
SN 1992A. FigurefTT]shows that within the luminosity region 
popu lated by SNe la which obey the Phillips relation the 
E08I data do not c onfirm the linear correlation suggested in 
Foley et al.l d2008l) . The range and scatter of Ruv measure- 
ments of the IEOSi sample appear to hav e more in common 
with the supernovae in the more recent IFolev et alJ (|2012l ) 
study which uses My so cannot be directly compared to 
these values. 

A larger sample of well-observed SNe may be needed to 
clarify this disagreement. In any case, Ruv is a flux ratio, 
not a broad-band colour measurement, so it is very sensi- 
tive to slight changes in the UV spectrum, which can be 
caused not only by different abundances in the ejecta, but 
also by small shifts in velocity. At 2770A, a change of 20A 
is equivalent to a change in velocity of « 2000 km s -1 . This 
is velocity change which could quite easily occur in normal 
SNe la. The measurements will also be strongly influence by 
the phase of the supernova as we can expect velocity shifts 
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Figure 11. i?-band magnitudes plotted against Ruv for the 
model s in this paper and objects investigated in iFolev etajj 
||2008| ) and [EOS! . No correction for light-curve shape or colour 
has been made to the data. 



as observed in the optical. The phase-dependence of the A2 
is shown in Figure [6] 

The best opportunity to use the UV for light curve cal- 
ibration is the strong trend between Lb and the (UV1 - 
b) colour as shown in Figure [5] The b filter corresponds to 
the optical region where supernovae are less susceptible to 
metallicity differences, and the features maintain a constant 
morphology despite differences in flux. Figure [5] shows that 
for the same luminosity the UV changes dramatically with 
metallicity, while Lb is basically constant. 

The UV1 region changes dramatically with both lumi- 
nosity and metallicity which drives the trend seen Figure 
[5] ft is possible that with optimisation, this colour could be 
used for lightcurve standardisation. (U-b) also shows a linear 
trend in the observed data as the U filter extends down as 
far as 3000 A, but at lower luminosity the models show the 
trend is not monotonic with metallicity. This would make 
any standardisation difficult. 



4.6 Future Extensions of This Study 

The observational data used in this paper are subject to 
large errors in both luminosity and Ruv- New observations 
with higher S/N in the UV would reduce the errors on Ruv- 
Observations of supernovae with lower luminosities would 
also enable a thorough comparison with the model data in 
these ranges. 

In terms of the modelling, we have used only one well- 
studied supernova as the basis of our models. While our 
spectra appear normal at a variety of luminosities and metal- 
licities, they are artificial in having been produced by ar- 
bitrary changes to the input parameters. With more well- 
studied supernovae with UV data, we would have more start- 
ing points for the creation of more accurate models covering 
a wider range of initial luminosities. This may then help to 
produce more accurate replicas of observed supernovae and 
reduce some of the disagreement seen between the models 
and the data. Velocity changes particularly affect the single 
feature and narrow-band measurements so sets of models 



which explore different parts of this parameter space would 
be especially helpful here. 

In the future it may also be possible to do a similar 
analysis for pre- maximum spectra where spectral signatures 
of the unburnt material from the progenitor are more easily 
identified in the optical. 



5 CONCLUSIONS 

We have used maxim um- light spectra for a sample of 9 SNe 
la obtained by E08 and compared them to a series of ID 
models produced by a MC radiative transfer code in order 
to test whether the variations observed the UV region of 
the spectrum can be explained by metallicity and luminos- 
ity changes alone. Our model spectra were initially based 
on one- well observed example SN 2005cf. We summarise our 
conclusions in the following points 

• Our models replicate well the range of broad-band UV 
colours seen in the observed spectra, showing that the dis- 
persion increases as one progresses to shorted wavelengths. 
The (UVl-b) colour appears to depend strongly on metal- 
licity, and it might have the potential for use in light-curve 
standardisation because it is less sensitive to very small 
changes in absorption features than narrow-band indica- 
tors(Sections I4TT1 and I4T5.2I) ; however this requires further 
observational data for testing at lower luminosities. 

• We observe that Ai and A2 can be interpreted as peaks 
caused by the reverse fluorescence of photons into the UV 
region of the spectrum, mainly from intermediate mass el- 
ements and chromium (Section I4.2[) . Both Ai and A2 move 
towards the blue with increasing metal content in the upper 
layers of the ejecta, but the effect is highly non-linear and 
may not be the best metallicity indicator. 

• We see that Ruv in our models has a low dependence 
on luminosity and the effect of the metal content is complex 
(Section 14.311. Hi gh-z data do not confirm the relation of 
IFolev et al.l (|2008l ). and so we do not support using this index 
for light curve standardization. 

• W e have used high-z Type la supernovae spectra from 
(jEOSl) and measured observed values of Ruv for 9 objects 
around maximum light. We see that the measured values 
agree with those in the models, although they are affected 
by large errors from noise in the spectra (Figure [7} . 

• We have performed an extensive search using the model 
spectra in order to investigate the possibility that these re- 
sults are driven by one or two elemental transitions, but the 
results have shown that the UV spectrum is far too com- 
plicated for this. In light of this we suggest that the use of 
Ruv is not appropriate for light-curve standardisation as 
the trends with metallicity and absolute magnitude are not 
linear (Section |4~4)| . 



ACKNOWLEDGEMENTS 

We acknowledge support from the Italian Space Agency 
(ASI) under contract ASI/INAF n. 1/009/10/0 and 
1/016/07/0. ESW and SH are grateful for the hospitality of 
MPA, Garching and the Weizmann Institute of Science, Re- 
hovot during stages of this work. Joint work by AG and MS 
is supported by the Weizmann-UK "making connections" 



12 E.S. Walker et al. 



programme. Joint work by AG and PAM is supported by 
a Weizmann Minerva grant. AG further acknowledges sup- 
port from the ISF, an ARCHES award, and the Lord Sieff 
of Brimpton fund. 



REFERENCES 

Altavilla G., Fiorentino G., Marconi M., et al., 2004, MN- 

RAS, 349, 1344 
Astier P., Guy J., Regnault N., et al., 2006, A&A, 447, 31 
Brown P. J., Roming P. W. A., Milne P., et al., 2010, ApJ, 

721, 1608 

Bufano F., Immler S., Turatto M., et al., 2009, ApJ, 700, 
1456 

Conley A., Howell D. A., Howes A., et al., 2006, AJ, 132, 
1707 

Conley A., Sullivan M., Hsiao E. Y., et al., 2008, ApJ, 681, 
482 

Cooke J., Ellis R. S., Sullivan M., et al., 2011, ApJ, 727, 
L35 

Ellis R. S., Sullivan M., Nugent P. E., et al., 2008, ApJ, 
674, 51 

Foley R. J., Filippenko A. V., Jha S. W., 2008, ApJ, 686, 
117 

Foley R. J., Filippenko A. V., Kessler R., et al., 2012, AJ, 
143, 113 

Garavini G., Nobili S., Taubenberger S., et al., 2007, A&A, 
471, 527 

Guy J., Astier P., Baumont S., et al., 2007, A&A, 466, 11 
Hamuy M., Trager S. C, Pinto P. A., et al., 2000, AJ, 120, 
1479 

Hoflich P., Wheeler J. C, Thielemann F. K., 1998, ApJ, 
495, 617 

Howell D. A., Sullivan M., Brown E. F., et al., 2009, ApJ, 
691, 661 

Howell D. A., Sullivan M., Conley A., et al., 2007, ApJ, 
667, L37 

Iwamoto K., Brachwitz F., Nomoto K., et al., 1999, ApJ 
SS, 125, 439 

Kessler R., Becker A. C, Cinabro D., et al., 2009, ApJ SS, 
185, 32 

Lentz E. J., Baron E., Branch D., et al., 2000, ApJ, 530, 
966 

Lucy L. B., 1999, A&A, 345, 211 

Mazzali P. A., 2000, A&A, 363, 705 

Mazzali P. A., Lucy L. B., 1993, A&A, 279, 447 

Mazzali P. A., Podsiadlowski P., 2006, MNRAS, 369, L19 

Milne P. A., Brown P. J., Roming P. W. A., et al., 2010, 

ApJ, 721, 1627 
Nomoto K., Thielemann F.-K., Yokoi K., 1984, ApJ, 286, 

644 

Oke J. B., Cohen J. C, Carr M., et al., 1995, PASP, 107, 
375 

Pastorello A., Taubenberger S., Elias-Rosa N., et al., 2007, 

MNRAS, 376, 1301 
Pauldrach A. W. A., Duschinger M., Mazzali P. A., et al., 

1996, A&A, 312, 525 
Perlmutter S., Aldering G., Goldhaber G., et al., 1999, ApJ, 

517, 565 

Perrett K., Balam D., Sullivan M., et al., 2010, AJ, 140, 
518 



Phillips M. M., 1993, ApJ, 413, L105 
Riess A. G., Filippenko A. V., Challis P., et al., 1998, AJ, 
116, 1009 

Riess A. G., Strolger L.-G., Casertano S., et al., 2007, ApJ, 
659, 98 

Sauer D. N., Hoffmann T. L., Pauldrach A. W. A., 2006, 
A&A, 459, 229 

Sauer D. N., Mazzali P. A., Blondin S., et al., 2008, MN- 
RAS, 391, 1605 

Stehle M., Mazzali P. A., Benetti S., et al., 2005, MNRAS, 
360, 1231 

Sullivan M., Conley A., Howell D. A., et al., 2010, MNRAS, 
406, 782 

Sullivan M., Guy J., Conley A., et al., 2011, ApJ, 737, 102 
Timmes F. X., Brown E. F., Truran J. W., 2003, ApJ, 590, 
L83 

Walker E. S., Hook I. M., Sullivan M., et al., 2011, MNRAS, 
410, 1262 

Wang X., Li W., Filippenko A. V., et al., 2009, ApJ, 697, 
380 

Wang X., Wang L., Filippenko A. V., et al., 2012, ApJ, 
749, 126 

Wood-Vasey W. M., Miknaitis G., Stubbs C. W., et al., 
2007, ApJ, 666, 694 



APPENDIX A: SINGLE ELEMENT 
SEQUENCES 

Figures [All - [A"5l show the effect that changing the amount 
of only 1 element has on the UV colours, positions of the 
UV emission features and the ratio Ruv- 
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Figure Al. The effect of changing just the amount of chromium in the outer ejecta. Upper - UV colours; Middle - The positions 
and A2; Lower: Rjjv 
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Figure A2. The effect of changing just the amount of stable iron in the outer ejecta. Upper — UV colours; Middle - The positions of Ai 
and A2; Lower: Ruv 
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Figure A3. The effect of changing just the amount of manganese in the outer ejecta. Upper - UV colours; Middle - The positions of 
Ai and A2; Lower: Ruv 



16 E.S. Walker et al. 



0.5 



-1.0 



i 5.00 
2.00 
1.00 

l 0.50 
0.20 

i 0.10 

i 0.05 




i 5.00 
2.00 
1.00 

l 0.50 
0.20 

i 0.10 

i 0.05 



UV2-b 




i 5.00 
2.00 
1.00 

l 0.50 
0.20 

i 0.10 

i 0.05 




8.6 8.8 



8.6 8.8 

logtivg 



8.6 8.8 
log(iyg 



l 5.00 
2.00 
1.00 

l 0.50 
0.20 

l 0.10 

l 0.05 




O Pre-max 
• Max 

c Post-max 



l 5.00 
2.00 
1.00 

l 0.50 
0.20 

l 0.10 

l 0.05 



> Pre-max 
i Max 

j Post-max 




' ■ 5.00 




2.00 




1.00 




- ■ 0.50 




0.20 




■ 0.10 




- m 0.05 










0.0 
8.4 



logflVLj 



Figure A4. The effect of changing just the amount of 56 Ni in the outer ejecta. Upper — UV colours; Middle - The positions of Ai and 
A2; Lower: Ruv 
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Figure A5. The effect of changing just the amount of titanium in the outer ejecta. Upper - UV colours; Middle - The positions of Ai 
and A2; Lower: R\jv 



